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1 Introduction 



Theories based on Supersymmetry (SUSY) are widely considered as the the- 
oretically most appealing extension of the Standard Model (SM). They predict the 
existence of scalar partners fi, fn to each SM chiral fermion, and spin-1/2 partners 
to the gauge bosons and to the scalar Higgs bosons. So far, the direct search for 
SUSY particles has not been successful. One can only set lower bounds of C(IOO) 
GeV on their masses p[. Contrary to the SM, two Higgs doublets are required in the 
Minimal Supersymmetric Standard Model (MSSM) resulting in five physical Higgs 
bosons . The direct search resulted in lower limits of about 90 GeV for the neutral 
Higgs bosons 

An alternative way to probe SUSY is to search for the virtual effects of the ad- 
ditional particles via precision observables. The most prominent role in this respect 
plays the p-parameter p. The leading radiative corrections to the p-parameter, Ap, 
constitute the leading process-independent corrections to many electroweak precision 
observables, such as the W boson mass, M^y, and the effective leptonic weak mixing 
angle, sin^^es- Within the MSSM the full one-loop corrections to and sin^ 
have been calculated already several years ago More recently also the leading 

two- loop corrections of 0{aas) to the quark and scalar quark loops for Ap and M]y 
have been obtained @,|]. Contrary to the SM case, these two- loop corrections turned 
out to increase the one-loop contributions, leading to an enhancement of the latter of 
up to 35% §. 

We summarize here the result for the leading two-loop corrections to Ap at 
0{Glmj) 0]. For a large SUSY scale, Msusy > Mz, the SUSY contributions 
decouple from physical observables. This has been verified with existing results at 
the one- loop ||lT| and at the two- loop level |]8|,p!0[|. Therefore, in the case of large 



Msusy the leading electroweak two-loop corrections in the MSSM are obtained in 
the limit where besides the SM particles only the two Higgs doublets needed in the 
MSSM are active. We derive the result for the 0{G'jpmf) |jlO| corrections in this case 
and provide a compact analytical formula for it, depending on the CP-odd Higgs 
boson mass. Ma, and the top quark mass, rrit. Furthermore, we present formulas for 
the limiting cases Ma ^ Trit (i.e. the SM limit) and Ma -C rrit. The numerical effect 
of the 0{G'jpmf) corrections is compared with the corresponding SM result |]12| and 
the gluon-exchange correction of 0{aas) in the MSSM. 



1 



2 Calculation of the 0{G%mf) corrections 

2.1 Ap and the Higgs sector 

The quantity Ap, 

Ap 



Ml 



parameterizes the leading universal corrections to the electroweak precision observ- 
ables induced by the mass splitting between fields in an isospin doublet 0. 'Lz^wiS^) 
denote the transverse parts of the unrenormalized Z and W boson self-energies at 
zero momentum transfer, respectively. The shifts induced by Ap in the prediction for 
the W boson mass, Miy, and the effective leptonic weak mixing angle, sin^ ^efr, are 
approximately given by 

Mw 



5M^ 



w 



r2 



r2 - q2 



Ap, (5 sin 6cs 



„2 2 
' r2 _ „2 
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Contrary to the SM, in the MSSM two Higgs doublets are required [Q. At the 
tree-level, the Higgs sector can be described in terms of two independent parameters 
(besides g and g'): the ratio of the two vacuum expectation values, tan/3 = V2/V1, 
and Ma, the mass of the CP-odd A boson. The diagonalization of the bilinear part 
of the Higgs potential, i.e. the Higgs mass matrices, is performed via orthogonal 
transformations with the angle a for the CP-even part and with the angle (3 for the 
CP-odd and the charged part. The mixing angle a is determined at lowest order 
through 



tan 2a = tan 2/? 



TT 



Ml' 



< a < 

2 



(3) 



One gets the following Higgs spectrum: 

2 neutral bosons, CP = +1 
1 neutral boson, CP = — 1 
2 charged bosons 
3 unphysical scalars : G'^,G^,G^. 

The tree- level masses, expressed through Mz, Mw and Ma, are given by 
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where the last two relations, which assign mass parameters to the unphysical scalars 
and G^, are to be understood in the Feynman gauge. 

2.2 Evaluation of the 0{G'^pm'f) contributions 

In order to calculate the 0{G'jpmf) corrections to Ap in the approximation that all 
superpartners are heavy so that their contribution decouples, the Feynman diagrams 
generically depicted in Fig. |l| have to be evaluated for the Z boson (V = Z) and the 
W boson {V = W) self-energy. We have taken into account all possible combinations 
of the t/b doublet and the full Higgs sector of the MSSM, see Sect. ^]T[ 

The two-loop diagrams shown in Fig. ^ have to be supplemented with the cor- 
responding one-loop diagrams with subloop renormalization, depicted generically in 
Fig. ^ The corresponding insertions for the fermion and Higgs mass counter terms 
are shown in Fig. |^. 

The amplitudes of all Feynman diagrams, shown in Figs. have been created 
with the program FeynArts2.2 [|1^], making use of a recently completed model file 
for the MSSM[[ The algebraic evaluation and reduction to scalar integrals has been 
performed with the program TwoCalc [H]. (Further details about the evaluations 



with FeynArts2.2 and TwoCalc can be found in Ref. IT^.) As a result we obtained 
the analytical expression for Ap depending on the one-loop functions Aq and Bq [|T6| 
and on the two-loop function T134 [p!4| , [l7| . For the further evaluation the analytical 
expressions for Aq, Bq and T134 have been inserted. In order to derive the leading 
contributions of 0{G'jpmf) we extracted a prefactor hf ~ G'jpmf. Its coefficient can be 
evaluated in the limit where Mw and Mz (and also nib) are set to zero. Furthermore 
we made use of the mass relations in the MSSM Higgs sector, see eq. (|^). In the limit 
Mw, Mz they reduce to 

ml = 
< = Ml 
= Ml 
ml = 

ml± = 0. (6) 

In the limit Mz — * the relation between the angles a and /3, see eq. (H), becomes 
very simple, a = j3 — tt/2, i.e. sin a = — cos (3 , cos a = sin (3 . The coefficient of the 
leading 0{G'j;^mf) term thus depends only on the top quark mass, m^, the CP-odd 

Higgs boson mass, M^, and tan (3 (or sp = tan j3/ \Jl + tan^ (3 ). 

We explicitly verified the UV-finiteness of our result. As a further consistency 
check of our method we also recalculated the SM result for the 0{Gpmf) corrections 

^ Only the non-SM like counter terms had to be added. 
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Figure 2: Generic Feynman diagrams for the vector boson self-energies with counter term 
insertion {V = {Z, W}, q = {t, b}, </>, x = {h, H, A, H^, G, 
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Figure 3: Generic Feynman diagrams for the counter term insertions 
{q = {t,b},<f> = {h,H,A,H^,G,G^}). 
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with arbitrary values of the Higgs boson mass, as given in Ref. [|T8|, and found perfect 
agreement. 



3 Analytical results 
3.1 The full result 



The analytical result obtained as described in Sect. ^]2| can conveniently be ex- 
pressed in terms of 



The corresponding two-loop contribution to Ap then reads: 
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with 

A = 3 - 13a + lla^ 
In the limit of large tan/? (i.e. (1 — s|) ^ 1) one obtains 
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Thus for large tan j3 the SM limit with Mf^ [jT2| is reached. 
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3.2 The expansion for large Ma 

The resuh for Ap^^^^g in eq. 
large values of Ma- 



can be expanded for small values of a, i.e. for 
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In the limit a — >^ one obtains 
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i.e. exactly the SM limit for Mf^ — is reached. This constitutes an important 
consistency check: in the limit a — the heavy Higgs bosons decouple from the the- 
ory. Thus only the lightest CP-even Higgs boson remains, which has in the OlCj^mf) 
approximation the mass rrih = 0, see eq. (P). This decoupling of the non-SM contri- 
butions in the limit where the new scale (i.e. in the present case Ma) is made large 
is explicitly seen here at the two-loop level. 

3.3 The expansion for small Ma 

The result for Ap^^lj^^^ in eq. (H) can also be expanded for large values of a, i.e. 
for small values of Ma (with a = 1/a): 
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In the limit a ^ or a ^ oo one obtains 
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(14) 



4 Numerical analysis 
4.1 The expansion formula 



We first analyze the validity of the two expansion formulas, eqs. (ITT]) and (0). In 



Fig. ||we show the result for ^i Hi^s^ defined by 

A „SUSY _ o '^F 4 .SUSY /I c-N 

^Pl,Higgs - 128 TT^ * "l.Higgs ' U^J 

as a function of 6 = MA/mt (= l/-\/a) fo^^ tan/3 = 3. The expansion for 6 <^ 1 is 
sufficiently accurate nearly up to 6 = 1. The other expansion gives accurate results 
for b >2. For larger tan (3 the expansion becomes better, enlarging the validity region 
for the large Ma expansion up to 6 > 1. 

4.2 Effects on precision observables 

In this section we analyze the numerical effect on the precision observables and 
sin^6'efr, see eq. (^, induced by the additional contribution to Ap. In Fig. |^ the size 
of the leading 0{a'^) MSSM corrections, eq. (H), is compared for tan/9 = 3,40 with 
the leading 0{a'^) contribution in the SM for Mf^ = 00, with the leading MSSM 
corrections arising from the t/b sector at 0{a) 0, and with the corresponding gluon- 
exchange contributions of 0{aas) |P (the 0{aas) gluino-exchange contributions [§], 
which go to zero for large rrig, have been omitted here). For illustration, the left 
plot (tan/? = 3) is shown function of Ma, which affects only the 0{a^) MSSM 
contributions, while the right plot (tan f3 = 40) is given as a function of the common 
SUSY mass scale in the scalar quark sector, Msusy, which affects only the 0{a) 
and 0{aas) MSSM contributions. We have furthermore chosen the case of "maximal 
mixing" in the scalar top sector, which is realized by setting the off-diagonal term in 
the t mass matrix, Xt, to Xt = 2Msusy and yields the maximal value for rrih for a 



given tan/5 (see Ref. [0 for details). In the right plot the case of no mixing, Xt = 0, 



is also shown. The mixing in the scalar bottom sector has been determined by using 
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Figure 4: The quality of the expansion formulas, eqs. (|Tl|) and (0), is shown as a function 
of b = MA/mt{= 1/^). 



a bottom quark mass of mb = 4.5 GeV, and by setting the trilinear couplings to 
Ab = At and the Higgs mixing parameter to /z = 200 GeV. The O(a^) contributions 
in the SM and the MSSM are negative and are for comparison shown with reversed 
sign. 

While for small values of Msusy the 0{aas) gluon-exchange contribution in the 
MSSM is much larger than the O(a^) contribution from eq. (|^) (note that in this 
region of parameter space the approximation of neglecting the scalar-quark contribu- 
tions in the O(a^) result is no longer valid), they are of approximately equal mag- 
nitude for Msusy ~ 500 GeV (this refers to both the no-mixing and the maximal- 
mixing case) and compensate each other as they enter with different sign. In this 
region the two-loop contributions are about one order of magnitude smaller than the 
0{a) MSSM contribution. For Msusy = 1000 GeV the leading MSSM 0{a^) contri- 
bution is about three times bigger than the 0{aas) gluon-exchange contribution in 
the MSSM. 

For small tan p (left plot of Fig. D and moderate Ma (Ma ~ 300 GeV) the 
new 0{a^) MSSM corrections are about two times larger than the leading (9(a^) 
contributions in the SM for M^ = 0. For large Ma the decoupling of the extra 
contributions in the MSSM takes place and the 0{a'^) MSSM correction approaches 
the value of the leading O(a^) contributions in the SM for Mf^ = 0, as indicated in 
eqs. (0), dH). For large tan/3 (right plot of Fig. D the 0{a^) MSSM correction and 
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Fi gure 5: The contribution of the leading O(a^) MSSM corrections, Higg's ' ^® shown 
as a function of Ma for tan /3 = 3 (left plot) and as a function of MsusY for tan /? = 40 
(right plot). In the left plot the case of maximal t mixing is shown, while the right plot 
displays both the no-mixing and the maximal-mixing case. Ap^^^^^" is compared with 
the leading 0{a^) SM contribution with Mff^ = and with the leading MSSM corrections 
originating from the t/b sector of 0{a) and 0{aas)- Both O(a^) contributions are negative 
and are for comparison shown with reversed sign. 
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Figure 6: The leading 0{a^) MSSM contribution to 6Mw (left plot) and (5sin^6'eff (right 
plot) is shown clS cl function of Ma for tan /3 = 3, 40. 
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the 0{a^) contribution in the SM for M|-^ = are indistinguishable in the plot, in 
accordance with eq. (|T0|). 

It is well known that the O(a^) SM result with = underestimates the result 



with realistic values of by about one order of magnitude |18]. One can expect 



a similar effect in the MSSM once higher order corrections to the Higgs boson sector 



are properly taken into account, which can enhance ruh up to rrih ^ 130 GeV [^], 
see Ref. 



In Fig. ^ the approximation formulas given in eq. (@) have been employed for 
determining the shift induced in Mw and sin^ O^s by the new O(a^) correction to 
Ap. In Fig. P the effect for both precision observables is shown as a fuction of Ma 
for tan/3 = 3,40. The effect on 6M-\^ varies between —1.5 MeV and —2 MeV for 
small tan/3 and is almost constant, 6M]y ^ —1.25 MeV, for tan/3 = 40. As above, 
the constant behavior can be explained by the analytical decoupling of tan/? when 
tan/3 3> 1, see eq. (10). The induced shift in sin^ 6'efr lies at or below 1 x 10~^ and 
shows the same qualitative tan/? dependence as SMyy. 



5 Conclusions 

We have calculated the leading 0{G'j^mf) corrections to Ap in the MSSM in the 
limit of heavy squarks. Short analytical formulas have been obtained for the full 
result as well as for the cases Ma ^ rrit and Ma rrit. As a consistency check we 
verified that from the MSSM result the corresponding SM result can be obtained in 
the decoupling limit (i.e. Ma oo). 

Numerically we compared the effect of the new contribution with the leading 0{a'^) 
SM contribution with Mf^ = and with the leading MSSM corrections originating 
from the t/b sector of 0{a) and 0{aas)- The numerical effect of the new contribution 
exceeds the one of the leading QCD correction of 0{aas) in the scalar quark sector 
for MsusY ^ 500 GeV. It is always larger than the leading C(«^) SM contribution 
with Mf^ = 0, reaching approximately twice its value for small tan /3 and moderate 
Ma. 

The numerical effect of the new contribution on the precision observables Mw 
and sin^ 6'efr is reletively small, up to —2 MeV for Mw and +1 x 10~^ for sin^6'efr. 
It should be noted, however, that the (9(q!^) SM result with Mf^ = 0, to which 
the new result corresponds, underestimates the result with realistic values of Mf^ by 
about one order of magnitude. A similar behavior can also be expected for the MSSM 
corrections. An extension of our present result to the case of non-zero values of the 
lightest CP-even Higgs boson mass will be undertaken in a forthcoming publication. 
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